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I. INTRODUCTION 


The purpose of this document is twofold. First, the document pre- 
sents a summary of the recent inner radiation zone electron data avail- 
able from the National Space Science Data Center. Second, the document 
presents an intercomparison of these data, describing the morphology of 
the inner radiation zone, and compares the data with the current genera- 
tion of inner zone trapped electron models. 

An analytic representation of the inner zone equatorial pitch angle 
distribution is presented. This model is based upon data from eight sat- 
ellites and is used to reduce all data to the form of equatorial flux. 

In general, in previous data comparisons, different pitch angle distribu- 
tions were used for data reduction; an important feature of the present 
analysis is the use of a common distribution. The data are compared with 
the inner zone models AE-6 ( Teague et al . , 1976) and AE-5 1967 ( Teague 
and Vette , 1972). For those L values for which storm-time variations are 
large, the data are event integrated using a simple storm model. 

The conclusion from this analysis is that, with few exceptions, the 
data are observed to be in good agreement. Although it is evident that 
no Starfish-free high-energy electron measurements are available from the 
inner portion of the inner radiation zone, the general conclusion is that 
the AE-6 model provides a good description of the present solar maximum 
environment. 


II. MORPHOLOGY OF THE INNER RADIATION ZONE 

Data from the following satellites were used in the present analysis: 
0V3-3, 0V1-13, 0V1-19, 0G0 5, 0G0 1, 0G0 3, 1963-083C, and 0S0 4. The data 
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are summarized in Table 1 and the time and energy coverage are shown in 
Figure 1, In this figure the energy interval of a channel is marked by 
a vertical bar and the nominal energy of a threshold channel is indicated, 
by a single upward arrow. The time span of a data channel is represented 

by the length of the horizontal bar together with the spacecraft designa- 

/ 

tion. In many cases the time coverage varies with the L value and the 
horizontal bars indicate that data were obtained at some L values. For 
the OV3-3 and 0V1-19 satellites, the energy interval of some channels over- 
laps. The overlap regions are indicated by cross hatching. In general, 
data were sorted into 0.1 L bins for L > 1.4 and 0.05 L bins for 1 < 1.4. 
Processing of the individual data sets is discussed in Appendix A. 

A. The Equatorial Pitch Angle Distribution Model 

In order to reduce the data in an efficient and consistent manner, an 
analytic equatorial pitch angle distribution was developed. The following 
form was adopted for fluxes within approximately two orders of magnitude 
of the equatorial flux: 

log 10 J = C for 90° > a 0 s a h (1) 

and ' = a 0 + aix + azx 2 + a 3 x 3 + a^x 1 * for 

a h > a 0 > a m 

where x = (90 - a 0 ) 2 

These equations define a pitch angle a h above which the flux is indepen- 
dent of equatorial pitch angle, a 0 . For equatorial pitch angles in the 
range a h to a m , a polynominal is used. In the region of the atmospheric 
cutoff for L > 1.4, a 0 < a m> a high order polynomial of the form 

log 10 J = a + b a 0 - 18 (2) 
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is used for which an atmospheric cutoff pitch angle, a c , is defined as the 
angle at which the flux is three orders of magnitude below J(a m ), and then 

a = log 10 J (a m ) - 3Y/(Y-1) 
and b = 3 a c 18 /(Y-l) 

where Y = (a c /a m ) 18 

For L < 1.4, the a 0 -18 expression is not used and equation (1) is extended 
to the range a h > a 0 > a c for which a c is defined as the pitch angle at 
which the flux is four orders of magnitude below the equatorial value or 
1 electron/ (cm 2 * s* sr*keV) , whichever angle is greater. 

The coefficients of the equatorial pitch angle distribution model are 
given in Table 2, and the family of curves given by equation (1) is compared 
for various L values in Figures 2 and 3 by normalizing at 90° equatorial 
pitch angle to a flux of 1 electron/ (cm 2 * s* sr*keV) . In addition, to facil- 
itate comparison, the curves in the upper panel in Figure 2 are normalized 
to a m (L = 2.4) using the expression 

a 0 (L) = a 0 (L) + a m (2.4) - a m (L) 

Norm 

and the curves in Figure are normalized to a m (L = 2.2) using a similar ex- 
pression. Note that neither a 0 nor C are given in Table 2, because these 
parameters represent the absolute flux at any L value and energy and do not 
relate to the shape of the equatorial pitch angle distribution. 

Over the range 1.4 < L < 2.2, two sets of coefficients for equation (1) 
were required to adequately describe the energy dependence of the equatorial 
pitch angle distribution. For higher and lower L values, insufficient data 
were available to observe any energy dependence. Similarly* for equation 
(2), energy dependence was not included in the coefficients. 

The coefficients of equation (1) were determined by weighted least 
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squares fit to data from the following satellites: 1963-038C, 0V3-3, 

0V1-13, 0V1-19, 0G0 5, 0G0 1, and 0G0 3. The most useful data sets in this 
respect were 0V1-13,,0G0 5, and, at low L values, 1963-038C. Where possible, 

the coefficients were determined from data obtained during magnetically quiet 

% 

periods. The analytic fits are compared with representative data in Figures 
4 through 40 for various L values and energies. With the exception of the 
normalized curves shown for 1963-038C and 0S0 4, the equatorial differential 
flux shown in these figures has units of electrons/ (cm 2 ‘S*sr«keV). For the 
curves shown in Figures 5, 6, and 9, the 0S0 4 data have been normalized to 
the 1963-038C data at the lowest B value of the OSO 4 observations. In Fig- 
ure 11 the 1963-038C data are normalized to Day 258, 1967. 

Most of the data shown in Figures 4 through 40 were used only in deter- 
mining the fit given by equation (1). The low pitch angle data, obtained 
near the atmospheric cutoff and used for determining the coefficients of 
equation (2), are shown in Figures 41 through 48 for L = 1.4 through 2.4. 
Normalized data are shown from the following satellites: OGO 5, 0G0 3, OSO 4, 

1963-038C, 0V3-3, and 0V1-13. Because of the general shortage of data in 
the region of the atmospheric cutoff, the data were normalized with respect 
to time and energy at a 0 = a m . Because the 0V1-13 satellite provided the 
most significant data set for a 0 > a m , the data were normalized to 500 keV 
and day 140 1968. The 1963-038C data are not shown explicitly, but separate 
fits were made to these data and the resulting curves are shown in Figures 
41 through 46. 

Weighted least squares fits were performed to the data and the fitted 
curves are shown in Figures 41 through 48. ' The smallest pitch angle for 
which these curves are shown corresponds to the model atmospheric cutoff 
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pitch angle, a c . Estimates of a c based upon the data are shown as hori- 
zontal bars. These cutoff angles are summarized in Figure 49, where the 
L dependence is shown for 1.2 < L < 2.4. Because the mirror point B val- 
ue is a more sensitive parameter in this region, this parameter is also 
shown and may be compared with the H mm curves for 50 km and 100 km, 
calculated using the Goddard Space Flight Center 12/66 magnetic field 
model updated to an epoch of January 1968. 

The data presented in Figures 4 through 48 indicated that, in gen- 
eral, the shape of the quiet time equatorial pitch angle distribution in 
the inner zone is described to better than a factor of two accuracy for 
all pitch angles by the model given in equations (1) and (2) and Table 2. 
The exceptions are at low L values (L < 1.35), high L values (L > 2.2), 
and at pitch angles within approximately 3° of a c for which data coverage 
is poor and considerable scatter is evident. In these regions the fit is 
accurate to no better than a factor of three. 

As noted previously, the fit is largely based upon quiet time data. 

Using 0G0 5 data. West et at. (1973) observed considerable variation in 
the shape of the outer zone pitch angle distribution correlated with mag- 
netic activity. In the inner zone, magnetic activity has a significant 
effect on flux levels for L > 1.7. Unfortunately, with the present data 
base, the satellites providing good pitch angle coverage at these L 
values (viz. 0V1-13 and 0G0 5) did not have good coverage during magnetic 
storms. However, a limited amount of data is available, and quiet time 
and storm time pitch angle distributions are compared in Figures 50 and 
51 at L = 2.0 and 2.4, respectively. These data indicate that within 
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the limited pitch angle coverage, no time dependence of the pitch angle 
distribution is observed, therefore the fit described in this section is 
applied to storm time and quiet time data. 

B. Data Intercomparison and Morphology 

The fit presented in Section II. A. was used to reduce the data to 
equatorial flux for the purposes of intercomparing data and examining 

t 

temporal variations. Flux vs time plots are presented in Figures 52 
through 118 for the following satellites: 0V3-3, 0V1-13, 0V1-19, 0G0 5, 

0G0 1, and 0G0 3. 0S0 4 equatorial flux radial profiles are given m 

Figure 119. The problems involved with processing the individual data 
sets using the model are discussed in Appendix A. 

The data are compared and discussed under the two general headings 
of long and short term temporal variations. The former category includes 
such effects as the decay of particles injected by the Starfish detona- 

f 

tion in 1962 and long term flux changes correlated with solar cycle vari- 
ations. The latter category refers exclusively to magnetic storm effects. 
These effects are rarely greater than a factor of two in the flux for 
L < 1.7. For instance, the largest magnetic storm observed in the data 
used in the present analysis occurred in May 1967 (peak daily mean Dst = 
-243) and, in Figure 112 at L = 1.6, it can be seen that the 0G0 3 data 
show flux variations of less than a factor of two at any energy. 

Quiet time inner zone spectra are presented in Figures 120 through 
134 for L values 1.2 through 2.4. Data are shown from time intervals for 
which data are relatively abundant and for which several data sets are 
available. The data shown are grouped into four general epochs (see Fig- 
ure 1) : 
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Oct. - Nov.' 1964 


OGO 1 


Pre-Sept. 1966 storm OV3-3, OGO 3 
Pre-June 1968 storm OV1-13, 0S0 4, OGO 5 
Jan. 1970 OV1-19 

The general agreement between the OGO 3 and the 0V3-3 quiet time 
data is good, particularly for the data taken near 0.5 MeV. The only 
discrepancy is evident over the L range 1.4 to 1.8 for which the high 
energy (1.33 MeV) .OGO 3 data are between a factor of two and five higher 
than the comparable OV3-3 data. A possible explanation may be that the 
proton correction given by Vampola (1975) for the high energy 0V3-3 chan- 
nels may be excessive. 

The general agreement among the quiet time data from OGO 5, 0V1-13, 
and 0S0 4 is good. Hie exceptions are the two lower energy channels of 
the 0V1-13 detector, which appear to be significantly lower than the 
comparable data for L > 1.3. In addition, at a few L values, the OSO 4 
data are significantly lower (e.g., L = 1.5) or higher (e.g., L = 1.9, 

E = 397 keV) than the OGO 5 data. However, the general good agreement 
of the OSO 4 data with other data may be taken as a measure of the ac- 
curacy of the pitch angle distribution over a wide range of equatorial 
pitch angles because the actual OSO 4 fluxes are in excess of two orders 
of magnitude below the equatorial perpendicular fluxes shown in Figures 
121 through 129. 

The general quiet time morphology of the inner zone over the inter- 
val from near solar minimum in 1964 to near solar maximum in 1969 is 
characterized in Figures 135 through 137. Three L values are shown cor- 
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responding to the inner edge of the zone (L = 1.4), the peak of the zone 
(L = 1.8), and the slot region (L = 2.4). Data were interpolated and 
normalized to three energies using the quiet time spectra and averaged 
where appropriate at any given time: 80 keV (0G0 1 and 3 at ~84 keV and 

0G0 5 at 79 keV), 500 keV (0G0 1 and 3 at 491 keV, 0G0 5 at 479 keV, 

0V3-3 at 475 keV, 0V1-13 at 500 keV, and 0V1-19 at 444 keV and 537 keV) , 
arid 800 keV (0G0 5 at 822 keV, 0V3-3 at 712 keV and 957 keV, and 0V1-13 
at 670 keV). 

A salient feature of the data shown in Figures 135 and 136 is the 
significant flux decrease observed at the higher energies over the time 
period 1964 to 1966. This is attributable to the decay of the artificial 
particles injected by the Starfish explosion in July 1962. The extent of 
the influence of these particles on the inner radiation zone electron 
flux was discussed by Teague and Stassinopoulos (1972) . They suggest 
that the longer lived artificial particles are at approximately 1.6 MeV 
at L = 1.3 and that the fluxes of these particles reach the natural 
background levels approximately 8 years after Starfish (i.e., beyond the 
time coverage of the present data) . With the exception of 0G0 1 data 
none of the data presented in this report are influenced by artificial 
particles for L s 1.7 or at any L value for E < 300 keV. To a signifi- 
cant extent, however, 0V3-3 and 0G0 3 data are influenced outside this 
range. The OGO 5, 0V1-13, and OSO 4 data appear to be almost completely 
free of the effects of artificial particle with the exception of the high 
energy channels (E > 822 keV) from OGO 5 at low L values (L = 1.5). The 
separation of the data into three distinct groups indicating various 
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phases of the decay of the Starfish particles can be clearly seen in Fig- 
ure 122. The relatively abundant artificial particles are at the higher 
energies such as those shown in Figure 138 for L = 1.5. The data in this 
figure were normalized to 1.5 MeV, and the 0V3-3 and 0G0 3 data were av- 
eraged. The AE-2 model is taken from Vette et at. (1966). The STARAD data 
from late 1962 were taken from West (1976) . The Starfish model of Teague 
and Stassinopoulos (1972) shows excellent agreement with the data at all 
epochs between 1962 and 1968. 

An additional feature of Figures 135 through 137 is the flux varia- 
tions at energies, L values, and times away from those .that are Starfish- 
influenced. The. variations are dramatized by the guidelines and appear 
to show a quiet time flux increase from solar minimum to solar maximum. 

The largest variation is a factor of four, although in some cases the flux 
increase is only slightly larger than the error bars. The effect is more 
pronounced at the peak of the inner zone than in the slot region. Obser- 
vations over a wider Starfish-free time period are required before this 
effect can be convincingly associated with solar cycle variations. In 
fact, in the region of the inner zone peak where the effect is largest, 
the distinction between quiet time and storm time variations becomes dif- 
ficult to make and is perhaps arbitrary. For instance, in Figure 57 at 
L = 1.8, it is probably not appropriate to identify the flux level at 
712 keV before the May 1967 storm as quiet because the effects of the 
earlier September 1966 storm may still be decaying. On the other hand, 
at L = 2.0 in Figure 59, the decay process from the September 1966 storm 
is clearly complete before the May 1967 storm occurs and the prestorm 
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flux for the later storm is significantly higher than the earlier storm.' 

i 

An additional effect is apparent in Figure 92, in which the 0G0 5 data 
show at all energies a significant flux decrease leading up to the October 
1968 storm. Although there is no significant enhancement of the flux a- 
round the time of the June 1968 storm, a reduction of the flux occurs 
following this storm. The conclusion is that it is not possible to sep- 
arate completely the inner zone flux into two independent components of 
storm time and quiet time flux. 

Seven magnetic storms were observed during the time period for which 
data are presented in this report. The timing of these storms in rela- 

J 

tion to the data is shown in Figures 1, 136, and 137, and the relative 
magnitudes of the storms as indicated by the Dst index are shown in 
Table 3. Figure 1 shows that the two large events of September 1966 and 
May 1967 and the small event of June 1968 were observed by two satellites. 
However, it should be noted that little data are available from the 0V3-3 
satellite for the May 1967 storm (e.g.. Figure 59) and from the 0V1-13 
satellite for the June 1968 storm. The peak storm spectra for those three 
storms are compared in Figure 139 for L = 2.4. An approximate quiet time 
spectra is shown for comparison. With the exception of the 300-keV and 
475-keV 0V3-3 data, plotted as different symbols, the various comparable 
data sets agree very well. A similar situation is evident at L = 2.0 
shown in Figure 140. It is clear from Figures 139 and. 140 that, relative 
to the quiet time spectra, storm flux enhancements are more pronounced 
at higher energies. Also, at the lower L value, the enhancement of the 
high energy spectra is significantly reduced compared to the higher L 
value, whereas at lower energies the enhancement remains relatively un- 
changed. 
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The general characteristic of magnetic storm effects in the inner 
zone is that they become larger at any given L value as the Dst index 
associated with the magnetic field changes becomes larger. The storm 
effects decrease rapidly for the lower inner zone L values. Significant 
storm time flux variations are not observed below L = 1.7, except for 
large events such as the May 1967 storm. These characteristics are sum- 
marized in Figure 141, which compares the quiet time and peak storm time 
radial profiles at 500 keV for the seven storms observed during the time 
period studied in this report. Comparing this figure with Table 3, it 
can be seen that the peak fluxes at any L value order with the peak 
hourly mean Dst index. That is, if an inner zone magnetic storm enhance- 
ment occurs, its magnitude correlates with the Dst index. However, it 
is not always true that if a Dst event occurs, an associated inner zone 
enhancement will occur. A further observation may be made concerning 
the rise time of the storm flux to the peak value. At high L values in 
the inner zone, the rise is rapid and occurs shortly following the peak 
Dst change. For instance, in Figure 63 the peak flux occurs on day 249, 
which is 2 days after the peak Dst change. At lower L values the rise 
becomes more gradual. In Figures 59 and 57 at L = 2.0 and 1.8, respec- 
tively, the peak fluxes occur on Days 265 and 280. At all of these L 
values a rapid flux change does occur on Day 249, but the magnitude of 
this change becomes a smaller percentage of the peak flux toward the 

lower L values. 

1 

III. COMPARISON OF THE CURRENT ELECTRON MODELS WITH DATA 
The inner zone electron models presently distributed by NSSDC are 
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AE-6 ( Teague et at,., 1976) for solar maximum and AE-5 1975 ( Teague and 
Vette, 1974) for solar minimum. Both models are intended to reflect the 
natural electron environment and do not include the effects of the arti- 
ficial particles injected by the Starfish event as the earlier model AE-5 
1967 ( Teague and Vette , 1972) did. Most of the data presented in this 
document were obtained around the epoch of 1967 and, at low L values, 
include the effects of the artificial particles. These data will be 
compared with the AE-5 1967 model, although this model is not presently 
distributed. 

The dependence of the model flux on B value is identical for the 
three models AE-6, AE-5 1967, and AE-5 1975. The model pitch angle 
distribution is compared with the present polynomial pitch angle dis- 
tribution in Figure 142 for L = 1.5- 1.8, and 2.4. In each case the 
curves correspond to 100-keV electrons, and the region near cutoff is not 
shown. The curves for L = 1.5 are typical of the L range 1.4 to 1.7. 

The model curve was based largely upon the 0G0 3 data ( Teague and Vette , 
1972) with the result that this curve xs flatter away from the cutoff 
region than the polynomial pitch angle distribution. This is discussed 
at length in Appendix A. 5. At these L values the model and the polyno- 
mial pitch angle distributions agree within a factor of two for near 
equatorial pitch angles (a 0 ^45° at L = 1.5, for instance) but larger 
differences are apparent' at lower pitch angles. The curves shown for L 
=1.8 are typical of the L range 1.7 to 2.0 and show agreement to better 
than a factor of two for all pitch angles away from the cutoff region. 

The curves shown, for L = 2.4 are typical of the high inner zone Lvalues. 
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The model curve was based largely upon the 0G0 3 data that provided poor 
pitch angle coverage near the equator at these L values. The curves 
shown in Figure 142 for L = 2.4 are normalized at 45°, and it is evident 
that the agreement is good for pitch angles below 70°, but that factor 
of three differences can be seen at the equator. For L values in the 
range 1.25 to 1.35, the AE-6/AE-5 model pitch angle distribution was 
determined largely from the 1963-038C data. The polynomial distributions 
are compared with these data in Figures 6, 8, and 11 where it can be 
seen that the agreement is good. As a result, the polynomial and the 
model pitch angle distributions for these L values also show good agree- 
ment (better than a factor of two) . The model and data cutoff pitch 
angles and B values are compared in Figure 49. The model cutoff B value 
or pitch angle corresponds to the value for which the 40-keV omnidirec- 
tional flux equals 1 electron/ (cm 2 *keV*s) . Although the data points given 
in Figure 49 correspond to a cutoff that is defined a little differently 
(see Section I I. A), the change of flux with B in this region is so rapid 
that these differences are insignificant, and the data and model may be 
compared. The cutoff given by the present polynomial occurs at a con- 
sistently higher B value than that given by the present electron models, 
although the data and model agree within the error bars for L > 2 and 
L < 1.4. It should be noted that these differences do not represent a 
major restriction of the model because they occur in regions where the 
flux levels are extremely low. ' 

In comparing the data and the present inner zone models it should 
be noted that the models are not intended to reflect the flux at a given 
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point in space at a given time, but rather are intended to give the user 
an estimate of the time-averaged environment at that point in space. The 
potential difference between these situations can be seen from Figure 141, 
in which AE-6 can be compared with peak storm flux data. At high inner 
zone L values, storm effects with characteristic times on the order of 
days result in flux variations of several orders of magnitude. The model 
AE-6 endeavors to represent an approximate average storm situation and 
may differ from the maximum and minimum observed flux at any given time 
by factors in excess of 10. As a result, the model is compared directly 

l 

with data only for those energies and L values for which short term tem- 
poral variations are small. The models AE-5 1967 and AE-6 can be com- 
pared with these data in Figures 120 through 126. In addition, the 
models are compared with 1963-038C data for L £ 1.3 in Figure 143. The 
data in this figure were either measured at the equator or near enough 
to the equator for the observed flux to be pitch angle independent. The 
general agreement between the data and the model is good. As noted in 
Section II.B, the higher energy data for L < 1.7 are influenced by arti- 
ficial Starfish electrons. The AE-5 1967 model was primarily a fit to 
0G0 3 and 0V3-3 data taken over the time period 1966-67. In Figures 120 
through 126, therefore, the AE-5 1967 model should be compared with data 
from this epoch. In general, for those energies for which AE-5 1967 is 
significantly different from the AE-6 model, the former model is signif- 
icantly higher than the comparable data. At the time of the generation 

I 

of the AE-5 1967 model, a provisional OV3-3 data set was used. Subse- 
quently, the calibration of the OV3-3 spectrometer was changed, result- 
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ing in a significant -reduction in the apparent observed fluxes ( Teague 
and Vette, 1974). The high energy 0G0 5 and 0V1-13 data obtained in 
1968 are more appropriately compared with the AE-6 model because this 
model is intended to be free of artificial Starfish particles. However, 
the Starfish model of Teague and Stassinopoulos (1972) predicts that the 
1.53-MeV and 2.82-MeV 0G0 5 data contain a significant artificial com- 
ponent for L £ 1.5, and it is reasonable to expect that the natural flux 
levels are below those measured by 0G0 5 in early 1968. The AE-6 model 
reflects this qualitatively, although quantitative confirmation of AE-6 
for E > 1.5 MeV and L s 1.5 must await data taken at later epochs. 

At energies for which AE-5 1967 and AE-6 are not significantly dif- 
ferent, the general agreement between the models and the data is good. 
However, at some L values (e.g., L = 1.4, Figure 124), the model is lower 
than the data by a factor of two for energies E < 150 keV. This results 
from the differences between the pitch angle distributions used for AE-6 
and the present polynomial expressions. In Figure 142, at L = 1.5, it 
can be seen that the polynomial expression differs from the AE-6 model 
by close to a factor of two at 45°, and that normalization of data from 
this pitch angle would result in equatorial fluxes based upon the poly- 
nomial curve approximately a factor of two higher than those resulting 
from the AE-6 model distribution. 

As noted earlier, for those L values and energies for which the tem- 
poral variations because of storm activity are significant, direct com- 
parison of the model and data is not appropriate. Comparison of the model 
with data integrated over a given storm, however, would show how accu- 
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rately the model would predict short term average fluxes in the inner 
zone. The following simple storm model is assumed for the storm time 
flux j (t) : 

j(t) = [Aj, - (j„ o - j qi )] exp (-t/T) + (3) 

t As shown in Figure 144, this model postulates a sharp rise from a 
prestorm flux j„ o to a peak flux j p = j,,^ + Aj, followed by an exponential 
decay with characteristic depletion time T to a post storm flux j q ^ . This 
model is quite appropriate for the higher L values, such as L = 2.4 (shown 
in Figure 63), but for the lower L values, such as L = 1.9 (shown in Figure 
58), the sharp rise has become rounded and the simple model is only an ap- 
proximation. .Equation (3) may be integrated over time T E (the event in- 
tegration time) to give the event integrated flux as: 

J E (T e ) = T_[Aj, + (j q - j_ )] [1 - exp(-T E /x)] + j (4) 

•T* 0 » I 

l E 

Estimates of the parameters in this equation were determined for all 
storms and all data sets for L a 1.7, and the event integrated flux was 
obtained for T E = 30, 90, and 180 days. Estimates of the error in J E 
were obtained from the error equation: 


( 

(T E /T)e" TE/T ' 

2 a q 2 + 0 P 2 ) 

a E 2 = 0„ 2 + (J E - j ) 2 { 

1 - 

- 

1 I 

lx 2 

(1 - e- T * /T >_ 

(Aj, + j„ o -j q , ) 2 J 


where a q , a p , and are the errors in the post-storm flux j q ^ , the peak 
flux j p , and the depletion time t, respectively. 


Event-integrated radial profiles for T E = 90 days are shown in Figures 
145 and 146 for SOO-keV and 800-keV electrons, respectively. In some 
cases interpolation was performed to reduce the data to the desired 
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energy. Data from the April 1965 storm are not shown at 800 keV because 
little 0G0 1 data were available at 1330 keV. Comparing Figures 145 and 
146, it is interesting to note that the event integrated data have quite 
different shapes from the peak flux data as a result of the pronounced 
decrease of depletion time T with increasing L value (t ~ 90 d at L = 1.8 
and t = 9 d at L = 2.4 and 500 keV) . As may be expected, the AE-6 model 
is a considerably better estimator of the event integrated flux than of 
the peak flux. At 500 keV, AE-6 is between a factor of two and three 
lower than the event integrated flux for the largest observed storm (May 
1967) . At 800 keV, in Figure 146, the corresponding factor is three to 
four. Note that the September 1966 storm may perhaps be regarded as a 
typical inner zone storm and that the model AE-6 gives the event integrated 
flux for this storm to within the minimum model error' (factor of two) at 
both 500 keV and 800 keV. For L < 2.2, this is the case for all energies 
as may be seen from the spectrum shown in Figure 147. For L > 2.2, the 
AE-6 model is less successful at high energies, as may be seen in Figure 
148 for L = 2.4. For E > 800 keV, the model is typically a factor of 
three lower than the September 1966 event- integrated flux and an order of 
magnitude lower than the large May 1967 event. AE-6 must, therefore, be 
regarded as low by more than the minimum model error for L > 2.2 and E s 
800 keV. The observed differences, however, are consistent with the con- 
fidence codes for AE-6 given by Teague et at. (1976). The event in- 
tegrated flux (equation 4) is a function of the event integration time 
T E . Arbitrarily, data are presented for a 90-day event integration time, 
and while this may be a reasonable minimum mission duration, the event 
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integrated flux and thus the model's ability to estimate it may vary 
considerably with T E . At low L values and/or low energies for which 
either the peak flux is within a factor of two or three of the prestorm 
flux or t is long with respect, to T E , the event integration time has rel- 
atively little effect upon J E . On the other hand, for small t, the effect 
may be substantial. This is shown in Figure 149 where J E is shown as a 
function of T E for 712 keV at L = 1.7 and 2.4. Within the error bars 
given by equation (5), the event integrated flux for other values of T E 
and E > 400 keV (30 s T E s 180 days) may be estimated from the expression 
J E al/T E N where N varies linearly with L values from 1.0 at L = 2.4 and 
0.5 at L = 1.9. 


IV. CONCLUSIONS > ' 

In comparing the inner zone electron flux measurements made by the 
0G0 1, 0G0 3, OGO 5, 0V1-13, 0V1-19, 0V3-3, 0S0 4, and 1963-038C satel- 
lites, the general agreement was found to be excellent. A number of sys- 
tematic differences were observed resulting from the decay of the arti- 
ficial Starfish particles and the long term variation of the flux at quiet 
times that are possibly correlated with solar cycle variations. A 
number of differences were observed between comparable data sets that are 
considered to be indicative of possible problems, with certain measure- 
ments. These are summarized as follows: (1) an apparent saturation 

effect in the OV3-3 data gave spurious pitch angle distributions for 300- 
keV and 475-keV electrons for L s 2.1 and for higher energy electrons for 
L £ 1.7 (Appendix A.l). An effort to correct this by omitting pitch 
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angles over which this effect occurred resulted in low equatorial fluxes 
at 300 keV and 475 keV (Section II. B); (2) the 0G0 1 and 3 data pitch 
angle distributions for all energies and 1.4 s L < 1.7 are flatter than 
those given by other data sets. Despite this, the data showed self- 
consistency, and special 0G0 1 and 3 pitch angle distributions were 
obtained to determine equatorial fluxes (Appendix A. 5). These fluxes 
showed good agreement with other data sets. However, the 1.33 MeV 0G0 3 
data appeared to be a factor of two to five higher than the comparable 
0V3-3 data. It was suggested that this discrepancy may result from the 
proton correction to the latter data set (Section II and Appendix A.l); 
and (3) the two lower energy channels of the 0V1-13 spectrometer gave 
measurements significantly lower than comparable data (Section II. B) . 

The trapped electron models were compared directly with the data 
for L < 1.6. In general, the agreement is good and better than the model 
confidence codes would suggest. However ■, little Starfish-free high 
energy (E > 1.5-MeV) data are available for comparison with the AE-6 
model for L < 1.5. In this region confirmation of the model must await 
more such data. At higher L values the data were event integrated for 
the various storms and compared with the AE-6 model and, in general, the 
agreement was good. Three problem areas were identified on the basis of 
comparison with the data: (1) the model atmospheric cutoff B value is 

too low by approximately 0.005 gauss at all L values (Section III); (2) 
for E < 150 keV and low L values, AE-6 is low by a factor of 2 to 2-1/2 
(Section IV); and (3) for E 2 800 keV and L > 2.2 the model is low in 
relation to the integrated data (Section ill). 
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APPENDIX A. DATA PROCESSING AND REDUCTION 


The data used in this report are summarized in Figure 1 and Table 1. 

In this appendix the processing and reduction of the individual data are 
discussed and archival flux vs time data sets are presented. As noted 
in the Introduction, all the data are reduced using the same equatorial 
pitch angle distribution at a given L value and energy. As a consequence, 
meaningful intercomparisons of the data may be made. 

Al. 0V3-3 Data 

The OV3-3 data were obtained using the Aerospace Corporation elec- 
tron spectrometer. The principal investigator was Dr. A. L. Vampola. 

The experiment was described by Vampola (1969). These data provided good 
time coverage throughout the last half of 1966 but provided only a few 
points in 1967. Good pitch angle coverage was given for L < 2.0. In the 
present analysis, these are the only data that provide significant measure- 
ments for E > 822 keV. 

The 0V3-3 data were processed using the method documented by Vampola 
(1975) . In the inner radiation zone the data were contaminated by pro- 
tons. Correction for the proton contamination was achieved using the ex- 
pression 0.8 P/(l - cos a), where P is the proton monitor count rate and 
a is the local pitch angle. Following Vampola (1975), all data for 
a < 25° were edited out. 

Examples of the 0V3-3 equatorial pitch angle distributions are pre- 
sented in Figures 10, 15, 20, 27, 30, 33, and 36, and some normalized 
small pitch angle data are shown in Figures 47 and 48. In the course of 
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processing the data, a saturation-like effect was observed for the low 
energy channels. An example is shown in Figure 150. In general, the level 
of the saturation-like effect was time and L value independent. In the 
subsequent processing of the 0V3-3 data, equatorial pitch angles above the 
value a, , at which this effect became apparent, were excluded. Values of 
a, are given in Table 4. 

The OV3-3 data were reduced using the model equatorial pitch angle 
distribution given in Section II. A, and the resulting equatorial flux vs 
time plots are shown in Figures 52 through 63 for the L range 1.3 to 2.4. 
For L > 1.8, three plots are given per L value, each showing three chan- 
nels of data. At lower L values the highest three energy channels are 
fully contaminated by protons, therefore only two plots are shown. 

A2. 0V1-13 Data 

The 0V1-13 data were obtained by the Air Force Geophysical Labora- 
tory (AFGL) electron spectrometer. The principal investigator was Dr. 

L. Katz. The experiment was described by Paotini et al. (1968); recali- 
bration of the experiment was described by Rothwell and Moomey (1972) . 

Data were available for, at most, 14 separate days over a short period in 
1968 (60 days) , but the pitch angle coverage was excellent for all L values 
1.25 < L < 2.4. In this L range the electron data were contaminated by 
omnidirectional protons with E > 75 MeV (RothueZl and Moomey, 1972). In 
an effort to correct for this contamination it was assumed that the count 
rate registered by the 900-keV electron channel on quiet days (Days 137- 
159 1968) resulted entirely from protons. An example of the resulting 
proton B distribution is shown in Figure 151, and the radial profile of 
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the proton correction at the minimum B observed at each L value is shown 
in Figure 152. The proton B distribution was approximated by two linear 
equations giving the proton count rate, P, as: 

log P = a B + b for B < B- 
10 ll Cr 

and a B + b + B r , (a - a ) for B > B rr 

For values of B for which these expressions gave count rates less than 
6.5/s a proton count rate of 6.5/s was assumed. Values of the constants 
are given in Table 5 and an example is shown in Figure 151. Rotbwell and 
Moomey (1972) give the proton count rate for the spectrometer as P = A J p 
(> 50 MeV) where J P is the omnidirectional proton flux and A = 0.25 cm 2 . 
The AP-8 trapped proton environment ( Sawyer and Vette, 1976) may be used 
to estimate J p ; the results are shown in Figures 151 and 152. As a re- 
sult of the favorable comparison it may be assumed that the present pro- 
ton correction is adequate. Examples of the proton-corrected equatorial 
pitch angle distributions are shown in Figures 4, 7, 17, 19, 22, 23, 25, 
29, 34, and 40. Normalized low pitch angle data are shown in Figures 41 
through 48. 

The 0V1-13 data were reduced using the model equatorial pitch angle 
distribution and the resulting equatorial flux vs time plots are shown in 
Figures 64 through 78 for the L range 1.2 to 2.4. 

A3. 0V1-19 Data 

The 0V1-19 data were obtained using the Aerospace Corporation mag- 
netic electron spectrometers. The principal investigator was Dr. A. L. 
Vampola. The instrument and the geometric efficiency factors used in 
processing the data were described by Vampola (1971) . The data coverage 
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spans the period from early 1969 to the beginning of 1970. 


In this report only some preliminaiy results of the data analysis 
are presented. Measurements of the fluxes within ±3° of the locally 
mirroring direction were selected so that the scatter of the data, caused 
by the finite apertures of the detectors, can be reduced. Because the 
correction of the pronounced proton contamination in the inner zone has 
not yet been determined, only data with energies below 537 keV and L s 1.8 
are presented. In this case the effect of contamination is expected to 
be minimal. Examples of this provisional pitch angle data are shown in 
Figures 26, 32, 38, and 39. When more information about the proton con- 
tamination and the angular response of the detectors becomes available, 
more 0V1-19 data can be used. 

The 0V1-19 data were reduced using the model equatorial pitch angle 
distribution, and the resulting equatorial flux vs time plots are shown in 
Figures 74 through 85 for the L range 1.8 to 2.4. 

A4. 0G0 5 Data 

The 0G0 5 data were obtained using the Lawrence Livermore Laboratory 
electron spectrometer. The principal investigator was Dr. H. I. West, Jr. 
and the experiment was described by West et at. (1969). These data pro- 
vided good pitch angle coverage for E < 822 keV and for L > 1.30. The 
time coverage was good at high innct zone L values but was successively re- 
duced for L < 1.8. Examples of the equatorial pitch angle distributions for 
0G0 5 are shown in Figures 8, 12, 13, 14, 16, 18, 21, 24, 31, and 37, and 
some normalized low pitch angle data are given in Figures 41 through 48. 
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The OGO 5 data were reduced using the model equatorial pitch angle 
distribution and the resulting flux vs time plots are shown in Figures 86 
through 98 for the L range 1.3 to 2.4. 

A5. OGO 1 and OGO 3 Data 

The OGO 1 and 3 data were obtained using the University of Minnesota 
electron Spectrometers. The principal investigator was Prof. J. R. Winckler. 
A reevaluation of the calibration of these spectrometers was performed by 
Teague (1970) and the resulting calibration constants appear in Table 2. 

The data from these satellites provide good time, L value, and pitch angle 
coverage, except for the data from late 1967, which have poor pitch angle 
coverage at high L values and little coverage at low L values. Because 
the intensity thresholds of the two higher energy channels were set at a 
high level, there is a small amount of high energy OGO 3 data, and the data 
provided by OGO 1 largely reflect artificial Starfish particles. 

Examples of the OGO 3 equatorial pitch angle distribution are shown 
in Figures 28 and 35. In comparing the shape of OGO 1 and 3 equatorial 
pitch angle distributions with those from other satellites it became ap- 
parent that the OGO 1 and 3 distributions for 1.4 ^ L s 1.7 showed a 
significantly slower cutoff toward the lower pitch angles. This effect 
is shown in Figure 153 in which a fit to the OGO 3 data is compared with 
the polynomial model given in Section II. A. The origin of the effect is 
unknown but it is observed to be systematic, however. As a result, specific 
OGO 1 and 3 equatorial pitch angle distributions were derived for the 
L range 1.4 to 1.7 for the purposes of reducing the data to equatorial 
flux. The form of the distribution is identical to the polynomial given 
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in Section 1 1. A, and the coefficients are given in Table 6. Note that the 
specific 0G0 1 and 3 distributions agree with the model distribution for 
fluxes within an order of magnitude of the equatorial flux. Normalized 
small pitch angle 0G0 3 data at high inner zone L values are shown in 
Figures 47 and 48. Reduced 0G0 1 and 3 equatorial flux vs time plots are 
shown in Figures 99 to 108 and Figures 109 to 118, respectively, for the 
L range 1.3 to 2.4. 

A6. 1963-038C Data 

The 1963-038C data were obtained by the Applied Physics Laboratory 
(APL) integral flux experiment. The principal investigator was Dr. C. 0. 

Bo strom, and the data were presented by Beall (1969). These data were used 
to a limited extent in the present study primarily for determining the 
equatorial pitch angle distribution for L < 1.4 because data from other 
satellites in this region were generally scattered. Normalized mirror- 
point distributions for these L values are shown in Figures 5, 6, 9, and 
11. At higher L values the 1963-038C data were used in the region of the 
atmospheric cutoff. Fits were made to the mirror-point distributions 
given by Beall (1969) and are shown in the equatorial pitch angle distri- 
bution plots given in Figures 41 through 46. Flux vs time curves are not 
presented for 1963-038C because these were given previously by Beall (1969) 
for specified B intervals. If desired, these curves may be reduced to the 
equator using the model given in Section II. A. 

A7. 0S0 4 Data 

The 0S0 4 data were obtained by the Lawrence Livermore Laboratory 
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experiment. The principal investigator was Dr. J. A. Waggoner. The ex- 
periment ,was described by Knox (1972) . Hie data were of limited value 
because the pitch angle coverage was poor and a large proton correction 
was required at low L value {Knox f 1972) . However, the data were useful 
in defining the atmospheric cutoff region of the equatorial pitch angle 
distribution. Normalized 0S0 4 data are shown in Figures 41 through 46. 
Flux vs time curves are not presented because only a few days of data were 
available. However, the equatorial pitch angle distribution model given 
in Section II was used to generate the 0S0 4 radial profiles in Figure 119. 
It should be noted that large factors, 0(10 2 ), were used to reduce these 

data to the equator. The error bars in these figures result from the time 

/ 

and pitch angle averaging and the proton correction. 
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Table 1. Summary of Data Used in the Present Analysis (concluded) 
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Table 2. Coefficients for the Analytic Equatorial Pitch-Angle Distribution 
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a 3 x 10 10 

a 4 x 10 11 * 
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Value 

(keV) 

a x x 10 4 

a 2 x 10 7 

(deg) ( 

(deg) 

(deg) 

1.2 

All 

110.0 

-402.5 

475.2 

-2456.0 

75.0 

- 

61.6 

1.25 

All 

26.9 

-122.6 

134.5 

-620.9 

77.8 

- 

55.7 

1.3 

All 

48.86 

-105:4 

9.56 

-265.1 

68.0 

- 

50.7 

1.35 

All 

-6.312 

-0.9535 

6.142 

-56.05 

69.0 

- 

46.9 
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-27.36 

20.51 
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48.0 

43.3 

1.4 

>500 
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-32.96 

22.83 

-61.52 

1.5 

<500 

12.75 

-24.56 

14.41 

-30.86 

70.0 

41.0 

37.2 

>500 

7.416 

-20.53 

12.85 

-28.83 


<500 

6.883 

-12.08 

6.055 

-12.06 

70.0 

36.0 

32.4 

/ 

I . o 

>500 

9.811 

-1-7.15 

8.415 

-15.27 


<350 

2.104 

-6.711 

3.548 

-6.686 

72.5, 

32.0 

29.0 

1.7 

>350 

-1.665 

-5.106 

3.33 

-6.5 

1.8 

<822 

-0.2274 

-4.143 

2.152 

-3.632 

75.0 

30.0 

26.1 

>822 

-0.94 

-5.627 

2.618 

-4.172 

i n 

<822 

-1.762 

-2.632 

1.397 

-2.238 


27.5 



>822 

-2.056 

-3.301 

1.722 

-2.702 

75.0 

23.7 

O A 

<900 

-2.098 

-1.777 

0.8971 

-1.372 

77.5 

24.0 

21.5 

Z. U 

>900 

-2.179 

-2.287 

1.088 

-1.573 


<900 

-2.242 

-1.52 

0.7384 

-1.054 

77.5 

23.0 

19.85 

2.1 

>900 

-2.599 

-1.687 

0.858 

-1.235 


<822 

-2.68 

-0.919 

0.498 

-0.6898 

80.0 

22.0 


2.2 

>822 

-3.42 

-0.3856 

0.3346 

-0.5526 

18.3 

i 

2.3 

All 

-2.519 

-1.448 

0.6945 

-0.8938 

85.0 

20.0 

16.9 

2.4 

All 

-3.298 

-0.9659 

0.5169 

-0.6551 

90.0 

18.5 

15.5 
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Table 3. Inner Zone Geomagnetic Storm ,Data 


Date of Storm 

Peak Hourly Mean 
Dst Index 

Peak Flux at L = 2.4 
and 500 keV* 

April 1965 

-185 

2.9 x 10** 

September 1966 

-229 

5.0 x 10 4 

May 1967 

-418 

1.9 x 10 s 

June 1968 

-94 

9.8 x 10 3 

October 1968 

-231 

1.5 x 10 s 

May 1969 

-133 

1.6 x 10 4 

September 1969 

-128 

7.5 x 10 3 


*electrons/cm 2 »s«sr»keV. 
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Table 4. 0V3-3 Saturation Equatorial Pitch Angles* 


L 

Value 


Energy 

(keV) 


300 

475 

712 

957 

1.3 

58 

62 

68 

N 

1.4 

52 

52 

62 

68 

1.5 

47 

47 

58 

N 

1.6 

40 

50 

58 

N 

1.7 

46 

46 

70 

N 

1.8 

47 

47 

N 

N 

1.9 

48 

50 

N 

N 

2.0 

50 

50 

N 

N 

2.1 

52 

55 

N 

N 


*In units of degrees. 

N: No saturation- like effect observed. 
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Table 5. 0V1-13 Proton Correction Coefficients 


L 

Value 

ai 

*2 

bi 

Bo 

1.2 

-30.03 

-61.6 

8.117 

-193 

1.25 

-18.35 

-43.85 

6.179 

-1885 

1.3 

-13.34 

-53.84 

5.339 

-202 

1.35 

. -12.98 

-123.42 

5.339 

-2194 

1.4 

-8.18 

-26.46 

4.549 

-1763 

1.5 

-8.9 

-45.48 

4.458 

-2 

1.6 

-10.97 

-32.33 

4.549 

-202 

1.7 

-11.31 

-28.97 

4.161 

-2205 

1.8 

-16.33 

-3.4 

3.978 

-169 

1.9 

-13.98 

-5.65 

3.456 

-1275 

2.0 

-27.34 

-7.05 

3.814 

-061 

2.1 

-40.71 

-2.49 

3.887 

-0625 

2.2 

-25.84 

- 

2.883 

- 

2.3 

-30.57 

-5.39 

2.436 

-045 

2.4 

-18.87 


1.87 

- 


logio P = aiB + bi for B < B Cr 

= a 2 B + bi + B Cr (ai - a. 2 ) for B s B c - 
Units: P counts/s 
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Normalized Electron Flux 
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Figure 6. Normalized Mirror-Point Distributions at L = 1.25 
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EQUATORIAL PITCH ANGLE DISTRIBUTION 
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EQUATORIAL PITCH ANGLE 
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EQUhTORIAL PITCH hNGLE distribution 



58 


EQUATORIAL PITCH ANGLE 
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Normalized Electron Flux 



Normalized Electron Flux 
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Figure 43. Normalized Low Pitch Angle Data at L = 1.6 
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Normalized Electron Flux 
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Normalized Electron Flux 
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Figure 45. Normalized Low Pitch Angle Data at L = 1.8 
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Figure 46. Normalized Low Pitch Angle Data at L = 1.9 
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Normalized 



Figure 47. Normalized Low Pitch Angle Data at L » 2.2 
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Normalized Electron Flux 
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Figure 50. Comparison of Storm and Quiet Time Pitch Angle Distributions at L = 2.0 



Equatorial Differential Flux [electrons/ (cm 2 s * sr keV)] 



Figure 51. Comparison of Storm and Quiet Time Pitch Angle Distributions at L = 2.4 
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Figure 52 (concluded)' 
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Figure 53 (concluded) 
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Figure 56 (concluded) 
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Figure 59 (continued) 
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Figure 61 f continued! 
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Figure 62 (continued) 
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FILIN' VERSUS TIME 
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Figure 63 (continued) 
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Figure 63 (concluded) 
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DRV OF VEAR 1966 





Equatorial Perpendicular Flux [electrons/ (cm 2 s • sr * keVH 



Figure 119. 0S0 4 Equatorial Flux Radial Profiles 
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Equatorial Perpendicular Flux [electrons/(cm 2 s-sr keV)] 
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Equatorial Perpendicular Differential Flux [electrons/(cm 2 • s • sr keV)] 
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Equatorial Perpendicular Differential Flux [electrons/(cm 2 s * sr keV)] 
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Figure 123. Quiet Time Inner Zone Spectra at L = 1.35 
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Figure 125. Quiet Time Inner Zone Spectra at L = 1.5 
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Figure 126. Quiet Time Inner Zone Spectra at L = 1.6 
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Equatorial Perpendicular Differential Flux [electrons/(cm 2 s • sr • keV)] 
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Equatorial Perpendicular Differential Flux [electrons/{cm 2 s sr keV)] 
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Equatorial Perpendicular Differential Flux [electrons/(cm 2 . s • sr • keV)] 
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Figure 129. Quiet Time Inner Zone Spectra at L = 1.9 
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Equatorial Perpendicular Differential Flux [electrons/(cm 2 • s • sr • keV)] 
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Figure 131. Quiet Time Inner Zone Spectra at L = 2.1 
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Equatorial Perpendicular Differential 
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Equatorial Perpendicular Flux [electrons/(cm 2 s • sr • keV)] 
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144, Storm Model 
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Figure 145. Event Integrated Radial Profiles E * 500 keV 
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Event Integrated Flux [electron$/(cm 2 • s - sr • keV)] 



Figure 147. Event Integrated Spectra at L = 1.8 
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Event Integrated Flux [electrons/! 


Event Integration Time = 90d 
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Figure 148. Event Integrated Spectra at L = 2.4 
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Event Integrated Flux [electrons/(cm 2 • s • sr • keV)] 



Figure 149. Comparison of Event Integrated Fluxes with Integration Time 
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Figure 151. 0V1-13 Proton B Distribution at 



Proton Count Rate (counts/s) 


t 



L Value 

Figure 152. Radial Profile of Proton Correction for 0V1-13 
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Figure 153. OGO 3 Pitch Angle Distribution 
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